How Can Macroscopic Models Reveal 
Self-Organization in Traffic Flow? 
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Abstract 

In this paper we propose a new modeling technique for ve- 
hicular traffic flow, designed for capturing at a macroscopic 
level some effects, due to the microscopic granularity of the 
flow of cars, which would be lost with a purely continu- 
ous approach. The starting point is a multiscale method 
for pedestrian modeling, recently introduced in [5] , in which 
measure-theoretic tools are used to manage the microscopic 
and the macroscopic scales under a unique framework. In 
the resulting coupled model the two scales coexist and share 
information, in the sense that the same system is simulta- 
neously described from both a discrete (microscopic) and a 
continuous (macroscopic) perspective. This way it is possi- 
ble to perform numerical simulations in which the single tra- 
jectories and the average density of the moving agents affect 
each other. Such a method is here revisited in order to deal 
with multi-population traffic flow on networks. For illustra- 
tive purposes, we focus on the simple case of the intersection 
of two roads. By exploiting one of the main features of the 
multiscale method, namely its dimension-independence, we 
treat one-dimensional roads and two-dimensional junctions 
in a natural way, without referring to classical network the- 
ory. Furthermore, thanks to the coupling between the micro- 
scopic and the macroscopic scales, we model the continuous 
flow of cars without losing the right amount of granular- 
ity, which characterizes the real physical system and triggers 
self-organization effects, such as, for example, the oscillatory 
patterns visible at jammed uncontrolled crossroads. 

1 Introduction 

Traffic flow along a road is ruled by complex dynamics, 
which involve mutual interactions among the vehicles. The 
effect of such interactions is a modiflcation of the traveling 
velocity of cars through successive slowing down and accel- 
eration maneuvers, resulting globally in self-organized flow 
patterns that can be clearly seen at large scales. The chal- 
lenging goal of mathematical models of road traffic should be 
devising methods able to reveal the spontaneous emergence 
of macroscopic self-organization from the microscopic inter- 
action dynamics. Indeed, from an engineering point of view, 
the macroscopic information is ultimately the most useful 
one for quantitative purposes, since it is generally more ro- 
bust, i.e., less prone to fluctuations, than the microscopic 
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one, and furthermore it refers to quantities, such as car den- 
sity and flux, directly measurable in practice. On the other 
hand, the actual physics of vehicular traffic pertains to the 
scale of single cars and may not be fully caught simply by av- 
eraged approaches, because the most interesting phenomena 
occur out of equilibrium. 

Mathematical models of vehicular traffic have been tra- 
ditionally focusing on just either scale of description. At 
the microscopic scale, see e.g., l4l [51 l9l 1131 1141 |20I . cars are 
modeled individually, usually as points, and their intercon- 
nected dynamics are formalized through systems of ordinary 
differential equations inspired by the classical framework of 
Newtonian mechanics: 
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where Xi = Xi{t) is the position of the j-th car at time t 
along a one-dimensional road and a; is a material model of 
its acceleration, which depends on the current position and 
speed of other neighboring cars (in most cases, of the head 
car only). Conversely, at the macroscopic scale, see e.g., 
[2] [S] |7| 1161 1171 118| , cars are assimilated to a continuum with 
density p obeying the conservation law: 



dp 8 , , 



(1) 



which expresses the conservation in time of the total number 
of cars, V being their average speed. If llTl is supplemented 
by a phenomenological relationship v = v{p) linking the 
speed to the local density {speed diagram), one obtains the 
so-called first order models. If, instead, a further equation 
is joined to (IlJ, expressing the balance of linear momentum: 

at ox 

a being now a material model for the average acceleration 
of cars, one obtains the so-called second order models. 

The aim of the present paper is to go beyond the just 
recalled dichotomy microscopic/macroscopic, taking a point 
of view centered around the concept of vehicular traffic as a 
complex system. Namely, a system in which the causes lead- 
ing to macroscopic observable outcomes have to be sought 
at smaller scales, particularly the microscopic scale of the in- 
teractions among vehicles. In more detail, our interest is in 
the possibility to reproduce, on the macroscopic flow, some 
self-organized effects typical of car-to-car interactions, such 
as, for instance, the traffic light effect. By this we mean the 
oscillatory pattern showing up when groups of cars from in- 
coming roads occupy alternately a junction, giving rise to a 
sequence of density wave packets along the outgoing roads. 
This usually occurs at jammed crossroads not regulated by 
any imposed red-green traffic light cycle |11| . 



In order to achieve our goal, we will evolve a multiscale 
method, based on the mathematical measure theory, re- 
cently proposed for pedestrian traffic modeling [5] . The key 
idea of the method is to describe the distribution in space 
of a population of moving agents by a measure, which can 
be understood as their mass or, alternatively, their proba- 
bility distribution if properly rescaled with respect to the 
total number of agents. According to the spatial structure 
of the measure, from such an abstract setting one can derive 
either a Lagrangian microscopic or an Eulerian macroscopic 
model: the former is obtained by a combination of Dirac 
delta's centered at the spatial positions of the agents, the 
latter by a density making the probability absolutely con- 
tinuous with respect to the Lebesgue measure. Furthermore, 
considering that a measure can feature both a singular and 
a regular part at the same time, one can also obtain mul- 
tiscale models in which the individual discrete/microscopic 
and continuous/macroscopic models coexist and automati- 
cally complement each other, without requiring a domain de- 
composition and transmission conditions between the scales 
at the interfaces as it happens instead in most multiscale 
models available in the literature (see e.g., | 10lll2iri5l ). 

The evolution of our method, that we present in this pa- 
per, is twofold: 

• we complement the multiscale coupling with a multi- 
dimensional coupling. Indeed, we model microscopic 
dynamics as genuinely two-dimensional, in order to re- 
tain the detail of possible lateral displacements of single 
cars, and the macroscopic dynamics as one-dimensional, 
for such a detail is actually not needed in a gross aver- 
age context. This also avoids the nontrivial imposition 
of boundary conditions along the lateral edges of the 
roads; 

• we perform a variable-in- space multiscale coupling by 
means of measures parameterized by the space variable, 
in order to account for the proper scale in the proper 
portion of the domain. 

With reference to the aforementioned case study of the 
traffic light effect, we demonstrate with our model that a 
self-organized flow cannot be reproduced by purely macro- 
scopic dynamics. The multiscale coupling, with the mi- 
croscopic scale possibly acting at the junction only, allows 
instead the macroscopic scale to inherit the right amount 
of granularity, which generates clearly visible density wave 
packets beyond the junction. 



2 The multiscale model 

In this section we derive our multiscale model, which is 
methodologically inspired by the one we developed in [5] 
for crowd dynamics. We first present it in full generality, 
considering a genuinely two-dimensional domain, then we 
specialize it to the coupling of one-dimensional macroscopic 
and two-dimensional microscopic dynamics. 

Let UP c M^, p = 1, 2, be the road taken by the p-th 
population of cars, agreeing that p = 1 labels the horizontal 
road and p = 2 the vertical one, see Fig.^ Then the domain 
of the problem is H := Q^ U Q? , particularly the junction is 
the subset 0} r\0? . Moreover, let Xj'*" £ QP be the position 




Figure 1: Two one-way roads with a junction (shaded) 

at time t of the i-th car of the p-th population, that we 
assume to be ruled by the following equation: 
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where, from left to right: 
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f2 — > R^ is the desired velocity of the p-th popula- 
tion, namely the one a car would choose if it were alone 
along the road. For simplicity, we assume that this ve- 
locity is constant in time, and, in particular, that v^ 
is rightward and v^^^, upward as indicated in Fig. nl 



1 

dcs 



• A'"'', (J = 1, 2, is the total number of cars in the g-th 
population; 

• K'"' : n X f2 — >■ K*^ is the interaction kernel, that is a 
mapping modeling how the presence of cars of the g-th 
population may induce the «-th car of the p-th popula- 
tion to modify its desired velocity. When q = p inter- 
actions are endogenous, otherwise they are exogenous. 

Model ([2]| is agent-based, because it takes invariably a 
Lagrangian viewpoint: cars are tracked one-by-one in space 
as time goes by. The way in which we now pass to a more 
general Eulerian model is by looking at car distributions in 
space rather than at cars themselves. For this, we formally 
consider the Xj'^^'s as random variables and we denote by 
^^ their probability distribution in M? . In practice, for all 
measurable subset E C Q, we have: 
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The fact that /ij is not labeled by the index i translates 
the simplifying assumption that cars are indistinguishable 
within the same population, so that /i^ is the common law 
of all of the XJ'^'s, i = l, ..., NP. 

Using stochastic Ito's calculus and ensemble averages, we 
get from ^ the following equation for the measure /ij' (see 
|19| for technical details): 



dt 



+ V.(/.fvf) = 
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with 
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KP9(x, y)<iA?(y|x), (4) 



where /ij (-Ix) is the probability law of the g-th population 
conditioned to the point x. 

In order to close model l|3]l-(|4l we may assume /ij(-|x) = 
/x^(-). In this case, it is interesting to note that the empirical 
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(<5 = Dirac delta) 



(5) 



is a solution to lIsl-Q if and only if the Xj'^'s solve |2]l. 
Hence model (111 is contained into the mathematical struc- 
ture Isl-Q. However the latter is much richer, indeed it 



also allows us to choose for /^^ a continuous distribution: 
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where C? is the Lebesgue measure in R-^ (viz. the area mea- 
sure) and pj = p^(x) : n — > [0, -foo) is the density distribu- 
tion function in space of cars of the p-th population. 

Describing the distribution of vehicles via measure llsl cor- 
responds to modeling cars as geometrical points with mass 
structure, thus emphasizing the intrinsic granularity of their 
distribution. Conversely, using measure (Iml corresponds to 
modeling cars as massless points distributed in space, thus 
focusing mainly on their continuous stream. If we plug sep- 
arately ([5]l and (Iml into (Isl, we get the following dual mi- 
croscopic/macroscopic representation (see [5] for technical 
details): 
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dt 



+ V ■ (p?vn = 
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where v^ is given by (Ul. The two dynamics can be prof- 
itably coupled, toward a multiscale model, by the following 
choice of the conditioned measure: 
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where d £ [0, 1] is a parameter fixing how much the coupling 
is biased toward the microscopic or the macroscopic descrip- 
tion. The advection velocity resulting by inserting (jsll into 
jijl is: 



vr(x) = v5^^(x) + E ^M ''^ E K'^'Cx, X^'«) 
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which practically realizes the coupling of the microscopic and 
macroscopic scale in ^\. Notice that 6 = corresponds to a 
genuinely macroscopic coupling, with the mass points Xj''' 
passively advected by the continuous flow. Conversely, (9 = 1 
corresponds to a genuinely microscopic coupling, with the 



density p^ passively advected by the velocity field generated 
by the points Xj'^. Finally, < 9 < 1 corresponds to a 
multiscale coupling, in which both the mass points Xj''' and 
the density p^ contribute actively to the resulting advection 
velocity. 

2.1 Car interactions 

The multiscale model consisting of Bl, n\, and Jsll still 
needs the interaction kernel K*"* to be specified. Inspired by 
the classical one-dimensional microscopic follow-the-leader 
model (see e.g., [l]): 



which integrated in time yields 
V 



x, = C- 



(Xi + l - Xi^ 



(q, 7>0), 



(r? := 0/7) 



for a suitable constant C (to be possibly regarded as the one- 
dimensional counterpart of the desired velocity) , we general- 
ize it to the present two-dimensional framework by setting: 



KP«(x, y) = -min 



^-f-^,A/^^|ls^,,(.)(y)r(x,y), 

(9) 



where, from left to right: 

• ??''' > is the interaction rate, namely a parameter 
measuring the intensity of the interaction between cars 
of the p-th and g-th population; 

• M^'^ > is a sensitivity threshold necessary for avoiding 
the singularity when y approaches x (remember that 
cars are modeled, in any case, as dimensionless points); 

• 1^; is, for a given set E, the indicator function of E 
(i.e., ls(x) = 1 if X G £, otherwise); 

• SRpq{x) C n is the interaction neighborhood, i.e., a 
subset of the ball centered at x with radius R^"^ > 0, 
where the presence of neighboring cars actually perturbs 
the car located in x; 

• f(x, y) := (y — x)/ |y — x| is the unit vector pointing in 
the direction of the relative position of the interacting 



We remark that the interaction neighborhood makes car 
interactions nonlocal in space; furthermore, whenever it 
does not coincide with the whole ball, it makes them also 
anisotropic. Coherently with the observation that cars are 
indeed anisotropic particles [6], we model the endogenous 
(p = q) interaction neighborhood as the right half-ball in 
the road Sl^ and the upper half-ball in the road f!^, for 
translating the fact that drivers normally look ahead only. 
Conversely, we model the exogenous (p ^ q) interaction 
neighborhood as the intersection of the right and bottom 
half-balls in the road Q^ and the upper and left half-balls in 
the road C^, for translating the fact that drivers look also in 
the direction of the incoming opposite fiow of cars, especially 
at the junction. 



2.2 Coupling 2D microscopic and ID 
macroscopic dynamics 

Model |4|-(l7|-|8]l-l|9l is fully two-dimensional. Nevertheless 
a two-dimensional description, while being possibly relevant 
at the microscopic scale for capturing the detail of lateral 
displacements of single cars, is actually less necessary at the 
macroscopic scale, where the interest is mainly in the prop- 
agation of density waves along the road length. For such a 
purpose, a one-dimensional description is clearly enough. 

Our measure-theoretic setting allows us to obtain eas- 
ily and rigorously such a multiscale multidimensional cou- 
pling. It suffices to assume that the density pj = Pj(x), 
X = (xi, X2), depends actually on one space variable only, 
particularly the one corresponding to the longitudinal direc- 
tion in either road: 



dpi 
8x2 

dxi 



Pt=Pt{^i) in f^^ 



p1 = p1(x2) in02 



The mathematical model remains formally the same, but for 
the fact that the advection velocity v^* in the equation for 
the density, cf. (ItI, has to be projected onto the longitudinal 
direction of either road. Specifically, if u^ denotes the unit 
vector in the longitudinal direction of QP (i.e., the horizontal 
unit vector 1 in Q^ and the vertical unit vector j in Cl'^), the 
2D microscopic- ID macroscopic model reads: 






(10) 



where ■ denotes the standard inner product in R^. 

Notice that, thanks to the fact that v^ is projected 
in one dimension after being computed as genuinely two- 
dimensional, multiscale interactions are automatically eval- 
uated with the correct dimensionality. The microscopic two- 
dimensional points Xj''' perceive a spatially two-dimensional 
density, simply constant in the transverse direction. On 
the other hand, the macroscopic one-dimensional density 
p^ is affected only by the longitudinal component of the 
interactions with the mass points, which are however two- 
dimensional. 



2.3 Variable multiscale coupling 

The parameter 9, which determines the multiscale coupling, 
is assumed to be constant in (jsll. This implies, in gen- 
eral, that the microscopic and macroscopic scales contribute 
jointly to the dynamics everywhere in the domain f2. How- 
ever, one may argue that their specific contributions are ac- 
tually relevant in different sub-domains. For instance, it 
might be important to retain the microscopic detail at the 
junction, where the actual granularity of the flow of cars 
plays a major role in shaping the exogenous interactions be- 
tween the crossing populations. Conversely, away from the 
junction, where only endogenous interactions occur, a gross 
continuous description may be appropriate. 



This amounts to converting 9 into a function ^(x) : Q — > 
[0, 1], so that the conditioned measure depends now explic- 
itly on the conditioning point x: 
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(■|x) = e(x)— 5^5^,,, + (1 - 9(x))p«£2, (11) 
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It is trivial to check that ( |11| is indeed a probability measure 
on R'^ for each x £ f2. 

The scale switch at the junction described above can be 
obtained, for instance, by means of the following function: 

flifxGV 

ex=lyx=<^ 12 

10 if X e s2 \ y, 

where y C R^ is a neighborhood of the junction Sl^ n Q^ 
(see Fig.H. 

3 Numerical simulations 

In this section we present some numerical simulations that 
show the potential of the multiscale multidimensional model 
described above. The numerical method used for solving 
l |10| , along with the multiscale algorithm, can be easily de- 
rived from those described in 5 . The main difference here 
is that the equations for the densities are one-dimensional, 
therefore the two-dimensional quantities computed by the 
microscopic part of the algorithm have to be projected to 
scalar quantities before they can be used by the macroscopic 
part. 

The two roads forming the computational domain are 
the stripes U^ = [0, 200 m] X [95 m, 105 m] and D,^ = 
[95 m, 105 m] X [0, 200 m]. Microscopic cars can freely move 
in the whole domain, while car densities are constrained to a 
numerical grid consisting of N nodes in each road. The space 
step is Ax = 200/ N. Desired velocities are v^^^ = (10m/s)i 
and v^^^ = (10m/s)j. The exponent 7 in the expression of 
the interaction kernel K^' is fixed to 7 = 1. 

Inflow boundary conditions are prescribed as follows. 

• Microscopic cars of either population enter the domain 
from xi = and X2 = 0, respectively, every Aij, = 0.9 s, 
randomly distributed along the width of the roads. 

• Car densities are imposed for all t > at the same 
boundaries: p^{xp = 0) = p;,, where the constant pt, is 
computed in such a way that the total macroscopic mass 
flowing into the domain up to the final time Tmax > of 
the simulation equals the total number of microscopic 
cars. Recalling that the p^'s are probability densities, 
this is written as: 



± max 
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Figure 2: Test 1: simulation with 8 = (genuinely maero- 
scopie dynamics). No particular effect is visible at the junc- 
tion 



Figure 3: Test 1: simulation with 9 = 0.7 (multiscale dy- 
namics). Traffic light effect beyond the junction is visible in 
both densities and microscopic cars 



At the end of either road (xp = 200 m) we assume that cars 
can freely exit with their desired velocity. At the macro- 
scopic scale, this is equivalent to imagine p^{xp) = for 
Xj, > 200 m and t > 0. 

Note that describing the roads as one-dimensional do- 
mains (rather than two-dimensional) allows one to avoid 
boundary conditions along the lateral edges, which are diffi- 
cult to treat analytically and can affect heavily the solution. 

Finally, the (variable) time step is chosen as: 

At = min{0.05s, AtcFLJ, 

where AtQpj^ is the time step imposed by the CFL condition 
at the macroscopic scale: 



AtcFL max 
p— 1, 2 

ipS [0,200 m] 



|(v?.u'')(xp)|=Ax. 



Remark 1. The above parameters, as well as those that will 
be specified in the subsequent numerical tests, are intended 
to be realistic but exploratory. In particular, they do not 
stem from any real measurement. 

Remark 2. As recalled in the Introduction, the simulated 
junction is not regulated by any give-way rule nor by red- 
green traffic light cycles. 

3.1 Test 1: Inducing the traffic light effect 
in the density 

In this test we aim at showing how granularity triggers self- 
organization, and, in particular, which effect it has on the 
macroscopic dynamics. Here we chose A'^ = 100, rj^^ = r)^^ = 
Im^/s, »?12 = ^21 ^ 35m2/s, H" = R^^ = 10 m, R^^ = 
i?2l = 20 m, M" = M22 = 15m/s, Afi2 = M^i = 50m/s. 

Figure [2] shows the outcome of the simulation with 6 = 0, 
i.e., with the macroscopic scale leading the dynamics of both 
point cars and densities. The latter flow quite uniformly and 
do not show any notable behavior at the junction. 

Figure [3] shows instead the outcome of the same simula- 
tion with 8 = 0.7. Here a clear traffic light effect appears 




Figure 4: Computational domain partitioned in two regions: 
the microscopic scale leads the dynamics in the black area V, 
the macroscopic scale in the remaining region (H^ Uil'^)\V 



beyond the junction. This is due to the fact that interac- 
tions among point cars are much different from those among 
densities: more granularity expedites the break of symme- 
try and triggers the typical alternate occupation of the space 
also at the macroscopic scale. 

3.2 Test 2: Space-dependent 6 

In this test we assume the scale parameter 8 to depend on 
X, as described in Sect. [273] In more detail, we choose 8 = 1 
inside a neighborhood V of the junction Cl^ D Cl^ centered 
at (100 m, 100 m) and ideally contained in a square of edge 
length 40 m, and 8 = elsewhere (see Fig. Wl. This is a 
natural choice since we expect that the granularity plays 
a role mainly at the junction. The other parameters are 
N = 200, r;il = r;22 = 7m2/s, r?" = r;^! = 27m2/s, i?" = 
r22 = 5ni, _Ri2 = _r21 = 10 m, M" = M22 = 20m/s, 
Ml2 = m22 = 40m/s. 

Figure [s] shows the outcome of the simulation at time 
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Figure 5: Test 2: 9 depends on x. The microscopic scale 
leads near the junction, the macroscopic scale elsewhere 




Figure 7: Test 2: Car density along the vertical road at 
time t = 23s with e(x) = (top), e(x) = 1 (middle), and 
6»(x) = lv(x) (bottom) 



car (not shown) density accumulates, due to strong interac- 
tions with Dirac delta's. Such an effect perturbs the density, 
especially where it should be ideally constant (e.g., before 
the junction), see also the same effect in Fig. 6.1(c) of [5]- 
This motivates the use of the microscopic scale only where 
it is really needed. To better quantify this qualitative dis- 
cussion, we plotted in Fig. [t] the one-dimensional density 
along the vertical road. Fully macroscopic dynamics do not 
show any appreciable effect either at the junction or beyond 
it. Conversely, fully microscopic dynamics generate irreg- 
ular waves beyond the junction, as well as a noisy density 
profile before it. Finally, mixed dynamics (fully microscopic 
near the junction, fully macroscopic elsewhere) produce a 
constant density before the junction and clear density wave 
packets beyond it. 



Figure 6: Test 2: 
everywhere 



S(x) = 1: the microscopic scale leads "RpfprpnrPS 



t = 23 s. Only the leading scale is plotted, namely point 
cars near the junction and densities elsewhere, but we re- 
call that both scales are present everywhere and evolve at 
all times, thus no transmission conditions are needed at 
the interfaces of V. Microscopic interactions trigger a self- 
organizing alternate passage beyond the junction, perfectly 
visible in terms of density waves, which persists also when 
the microscopic scale ceases its influence on the dynamics. 
This test shows that self-organization can be obtained by 
adding granularity only at the junction rather than from 
the very beginning like in Test 1 (cf. Fig.jsll. 

It is interesting to compare the previous result with the 
extreme cases 6(x) = and 6(x) = 1. The simulation with 
vanishing 8 does not show any particular effect at the junc- 
tion, analogously to the result depicted in Fig. [2] Conversely, 
with 6{x) = 1 (dynamics fully led by the microscopic scale) 
the trafSc light effect comes up at and beyond the junction, 
see Fig. |6] However, an annoying effect also appears before 
and, less visible, beyond the junction: behind every point 
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